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Abstract. We present a method of directly obtaining the parity of a Gaussian 
f^ ' state of light without recourse to photon-number counting. The scheme uses only a 

^pH, simple balanced homodyne technique, and intensity correlation. Thus interferometric 

y ' schemes utilizing coherent or squeezed light, and parity detection may be practically 

jrt ' implemented for an arbitrary photon flux. Specifically we investigate a two-mode, 

^ I squeezed- light, Mach-Zchnder interferometer and show how the parity of the output 

state may be obtained. We also show that the detection may be described independent 

of the parity operator, and that this "parity-by-proxy" measurement has the same 

signal as traditional parity. 
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1. Introduction 

Quantum metrology, in part, is the investigation of methods of accurate and efficient 
phase estimation. To this end, various interferometric setups have been proposed. They 
mainly differ in the light they use and - as a consequence of this - the detection scheme 
that is required to obtain the phase information. [1] 

Parity has been shown to be a very desirous method of detection in interferometry 
for a wide range of input states. For many input states parity does as well, or nearly as 
well, as state-specific detection schemes [2] . Furthermore, as has been reported recently, 
parity paired with a two-mode, squeezed-vacuum interferometer actually reaches below 
the Heisenberg limit, achieving the Cramer- Rao bound on phase sensitivity. [3]. 
Mathematically, parity detection is described by a simple, single-mode operator 



n = f-i 



\N 



where N is the photon number operator. Hence, parity is simply the evenness or oddness 
of the photon number in an output mode. 

A potential advantage of parity detection might is metrology in the presence of 
loss. Many of the states which show the greatest potential phase sensitivity also quickly 
degrade in lossy environments, limiting their usefulness. In a real-world application 
of super-sensitive remote sensing it may, or may not, be advantageous to use delicate 
quantum light - depending on the current conditions in the environment. One could 
imagine an adaptable device which could be used to send one of several different states 
out into the environment to perform measurement. Perhaps a NOON state [1] for low 
loss, an MhM' state [5] for intermediate loss, and a coherent state for high loss [6l [7] . 
The receiver for such a device would not need to be different for each kind of light, but 
could always use parity detection, adding some robustness and ease of implementation. 

Experimentally, parity is often measured, directly, by counting the number of 
photons in the light field. Unfortunately, number-resolving detectors are difficult 
to make and operate, and only work at low photon numbers [8]. Though there is 
significant effort and continual progress in the field of quantum optical metrology 
P [TOl [m [121 [ISl im CSl ESI E], it wouM be useful to have another way of determining 
parity. Moreover, photon number counting is more than is necessary to obtain the parity 
of a state. Alternatively a non-linear optical Kerr interferometer may be used to find the 
evenness or oddness of a light field [17]. But the feasibility of this scheme rests on the 
availability of materials with large non-linearities. In this paper we present a method of 
obtaining parity directly - without recourse to cumbersome photon number resolving 
detectors or Kerr materials. Our method is much simpler for states of light, which 
have Gaussian Wigner functions (a class which includes both squeezed and coherent 
light). We focus on the two- mode squeezed- vacuum MZI mentioned above and present 
a complete setup which allows implementation of our promising scheme. 
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2. Review of Two-Mode Squeezed Vacuum Fed Mach-Zehnder 
Interferometry 

In a recent paper by Anisiniov, et al. [3] , it was shown tliat a MZI fed witli a two- mode 
squeezed vacuum state, and utilizing parity detection, could reach below the Heisenberg 
limit on phase sensitivity in the limit of few photons. In the high-photon limit the 
sensitivity goes as A0 ^ l/fi, approximately the Heisenberg limit. For all photon- 
numbers the setup hits the Cramer-Rao bound. The setup is also super-resolving. A 
diagram of this kind of setup is provided in Fig. [TJ 




Figure 1. A diagrammatic representation of a two-mode squeezed vacuum fed MZI. 
An optical parametric amplifier is pumped with a strong coherent beam, emitting 
highly correlated photons into two vacuum modes. These modes then interact with a 
50:50 beam splitter, a phase shift, and another beam splitter. Parity detection may 
then be performed on either mode. This setup is both super-sensitive and super- 
resolving. 



3. Practical Measurement of Parity Using Balanced Homodyne 

The relationship between parity and the value of the Wigner function at the origin has 
been known for some time |18]. It is simply 



n 



TT 



W^(0,0). 



(1) 



We present an entire and concise derivation of this fact for the sake of clarity 
and completeness in the Appendix. Wigner functions of unknown quantum states 
are typically reconstructed after performing optical quantum state tomography [T9] . 
Homodyne measurements are used to measure the the expectation values of the creation 
and annihilation operators at different bias phases (each phase constituting a slice). The 
full Wigner function can then be built up from these images. This would be sufficient 
to discover the value at the origin and thus the parity, but it is overkill. We only need 
the value at one point. 

Let us consider a specific class of states, called Gaussian states. These states are so 
called because their Wigner functions are Gaussian in shape. This class includes a very 
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broad category of states, including the coherent states and, conveniently, the squeezed 
states. It has been shown [20] that Gaussian Wigner functions are given by 



1 _ u{a — ao) +u* {a — ao)* +t|q — Qo|' 

W{a,a*) = — e ^^^^^ 



where (a) = Oo, (a^^) — (d"'")^ = —2m, and {d^a+^) — {a'^){a) = r. So we need to find Oq, u, 
T, and their complex conjugates. Some examples: for coherent states u = 0, and r = 1/2, 
for two mode squeezed vacuum r = 2sinh^(r) + 1/2, and 2n = cosh(r) sinh(r), where 
r is the gain. Also, since squeezing the vacuum state does not change the expectation 
values of its quadratures Oo = 0. Hence, we are left with 

W{0,0)= / (2) 

tta/t^ -4 nP 



So now we need only find (d^d), which is simply the intensity of the light. We also 
need to find (d^^) or its complex conjugate, which is the expectation value of a non- 
Hermitian operator. Thus it can not be measured directly. However the technique of 
balanced homodyning may be employed to determine it indirectly. 

Our analysis is based on the fact that the Gaussian character of the input state 
does not change after it has gone through the interferometer, as the MZI is a linear 
device. 

To begin let's take the standard example of balanced homodyne measurement, 
originally developed by Yuen and Chan [21]. The problem that homodyne techniques 
address is how to measure an, usually inaccessible, quadrature of the light field. The 
unknown light beam (represented by the state vector {ip)), on which measurements are 
to be performed, is mixed on a beam splitter with a strong coherent beam of known 
intensity and phase (|/3)), called a local oscillator. Photo-detectors are then placed at 
the output ports. See Fig. |2l 

For the case of a 50-50 beam splitter we utilize the standard transformations 
c — )■ (d -|- ib)/\/2 and d ^ {b + id)/\/2. The intensity difference at the detectors in 
terms of the input operators is then given by 

ND = c^c-d^d=-i{d^b-b^d). (3) 

The expectation value of Eq. (|3]), taken at the input is 



mmowm = -m{{d^)e'' - {d)e- 



id 



where 6 is the known phase of the local oscillator. Let's call this expectation value Y{6). 
Now we see that we can measure either (d^) or (d) of the unknown light field by making 
measurements at two different phases of the local oscillator 
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1/3) 




Figure 2. An unknown light beam lip) is mixed on a beam splitter with a 
strong coherent beam of known intensity and phase |/3), called a local oscillator. 
Photodetectors are then placed at the output ports. These are connected to a post- 
processing unit which carries out a intensity differencing. Thus the quadratures of the 
unknown light may be measured. 
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Figure 3. The setup we will discuss: a strong pump beam is sent in from the left. 
Some of the beam is peeled off for later use as a local oscillator. The beam then goes 
through a frequency doubler before entering the TMSV MZI depicted in Fig. [TJ Thus 
the squeezed light and the reference beam arc of the same frequency, and are phase 
locked. On output, one mode is mixed on a beam splitter with the local oscillator. 
After this, two detectors make intensity measurements on the final outputs which are 
fed into a post processor, which controls the bias phase, and beam chopper placed in 
the local oscillator beam. 



W nloi 



Closely related to this is our problem of measuring the second order moments of 



Parity Detection in Quantum Optical Metrology Without Number Resolving Detectors 6 

squeezed light in order to obtain the Wigner function at the origin - and thus parity. 
We will need the second order moment (a^^). To do this take the setup in Fig. [31 A 
strong pump beam is sent in from the left. Some of the beam is peeled off for later use 
as a local oscillator. The beam then goes through a frequency doubler before pumping 
an OPA. Thus the squeezed light and the reference beam are of the same frequency, and 
are phase locked. The squeezed light then proceeds through a standard MZI, interacting 
with the phase to be measured. On output, one mode is mixed on a beam splitter with 
the local oscillator. After this, two detectors make intensity measurements on the final 
outputs which are fed into a post processor, which controls the bias phase, and beam 
chopper placed in the local oscillator beam, according to a prescription which will be 
discussed shortly. This scheme is designed to measure the second order moments of the 
lower output port of the MZI. 

First we guess the moment d'^dc^c, which represents the intensity-intensity 
correlations between the two detectors. When this is propagated back through the 
beam splitter it becomes 

d^dc^c ^ ^ {dfa} + b'df + P'd} + bH' 
Taking the expectation value, we define the function 

X{e,\P\)=4{d^dc^c) 



which constitutes a single measurement. It is (al ) that we wish to obtain. We can 
accomplish this by performing three X measurements at, 6' = 0, 6' = 7r/4, and |/3| = 
(i.e. when the beam is blocked), and arranging them according to the prescription 

(af) = ^ [^X{0, m) + X(vr/4, |/3|) - (z + 1)X(0, 0) - (z + 1)1/3^] (4) 

We can obtain {dldf) easily with (d^d) + {c^c) — |/3p. With this information and Eq. 
(|2]) we can obtain the parity of a TMSV fed MZI. 

But, does this detection scheme really produce the same signal as parity? To 
answer this question we can remove any mention of parity from the calculation and 
directly compute the signal of the detection protocol. In order to do this we employ 
an operator propagation technique. The operators at the output are related to the 
operators at the input by the matrix transformations in Eq ([5]). The first matrix (from 
right to left) represents the Bogoliubov transformation of the OPA, where fi = cosh(r), 
u = sinh(r), and r is the gain. The phase of the pump has been set to zero. The 
next four matrices represent the first beam splitter, the probe and control phases, the 
second beam splitter, and the homodyning beam splitter. We then can write down the 
the output operators, which we use during detection, in terms of the input operators. 
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where taking the expectation values is more tractable. Despite this a specially written 
computer code is still required to compute these expectation values. 
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A single measurement of X evaluates to 
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(5) 



X{9, 



1 

16 



11 + 16/3* + cos(20) - 16cosh(2r) 

+16/3^(cos(2^ - 0) sin(0) sinh(2r) 
-(cos(20) - 5) cosh(4r)] . 

Performing the three measurements and using Eq. (|4]) we obtain 



lt2\ 



) = e ^"^ cosh(r) sinh(r) sin(0). 

For the sake of simplicity we take the device to be lossless, making (n) = 
the intensity output of a TMSV fed MZI is independent of phase, that is 



sinh 



(6) 



r , as 



(outlaws flout) = (out|a|jan|out) = sinh^(r). 

There is a quick and intuitive proof of this somewhat surprising outcome. It is due to 
the fact that when a two-mode squeezed vacuum state, |.^tm), is incident on both ports 
of a 50:50 beam splitter, the output is the product of two single-mode squeezed vacuum 
states, |^c)|^d)- Thus, for beam splitter, and phase shift transformations: B, and P, 
the intensity of output mode A can be written as 



Ia = {^TulB^PlB^A^ABPcmTu) 

= {^c\{^d\C^C + D^D - ib^Ce'^ + iCWe-"^\ic)\iD) ■ 

Only the last two terms carry phase information and, since for squeezed vacuum states 
the expectation values of first order moments is zero, the phase dependence is eliminated. 
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Using Eq. ([2]) we can write the signal of the detection protocol, S, in terms of these 
operators as 



S 



Substituting our calculated expectation values, Eq. ([6]), and making the choice of bias 
phase (j) ^ (j) + 7r/2 this expression becomes 

s- 



l + n(n + 2)sin^(0) 

Which is identical to the signal for parity detection given in Ref. [3], where n is the 
total average number of photons exiting the OPA. Thus we conclude that our homodyne 
technique reproduces exactly parity detection. Though we have not computed explicitly 
the quantum noise for this setup, it is reasonable to assume that since the signal is 
precisely the same the noise will be equivalent. The calculated minimum detectable 
phase shift for the Anisimov, et al., setup is given in Ref. [3] as 

A0min - 



ln[n + 2) 
in the vicinity of = 0. 

4. Alternate Setup 

We would also like to present an alternate setup for obtaining the parity of a TMSV 
MZI. This setup is presented in Fig. HI 

Here, instead of performing parity by proxy on only one port we perform it on both 
simultaneously. This will allow noise, which is endemic to the device, to be compensated 
for. Furthermore this may provide some robustness to loss. Also, note that the beam 
choppers have been removed, this is because this device demonstrates an alternate way 
of eliminating the unwanted terms using only the controlled phase shift. Not using 
the beam chopper has the advantage of keeping a more consistent intensity on the 
photodetectors, which may allow for more sensitive devices to be used. The price of 
their removal however is that four measurements at four different settings of the bias 
phase are required (as opposed to only three in the previous setup). 

Let us take the bottom port, and redefine the X measurement as 

= {a{]a\,) + |/3| V^^(ag) + |/3| V^^^(a?,) + \P\\ 
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Figure 4. The alternate setup: a strong pump beam is sent in from the left. Two 
beam splitters peel of some of the pump beam for use later as local oscillators for 
each output port. On output, both modes are mixed on a beam splitter with local 
oscillators. After this four detectors make intensity measurements which are fed into 
post processors which control the bias phases placed in the local oscillator beams and 
computes parity in both ports. Also, note that the beam choppers have been removed. 



where we are now no longer concerned with the intensity of the local oscillator. We can 
obtain the desired moment by performing this measurement four times, at four different 
phases, according to the prescription 



t2, 



(alf) 



zX((0)+X((f -zX((f -X; -f 



4^|/3i 



Again this can be used to obtain parity, likewise for the upper port. 



5. Summary 

To conclude we have devised a method by which the very desirable parity measurement 
may be performed on Gaussian states. We have shown in detail how to set up 
a parity detector for the specific case of squeezed light using homodyning, thus 
realizing the possibility of practical sub-Heisenberg phase estimation in a Mach-Zehnder 
Interferometer. We also showed how, though the detection scheme mimics parity 
measurement, it may be considered conceptually independent and achieves exactly the 
same signal as parity. 

It is also useful to point out that the parity of Non-Gaussian fields may be obtained 
via Eq. ([T]), however in this case a measurement of the full Wigner distribution, using 
quantum state tomography, is required. 
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Appendix 

First, start with the characteristic function of the Wigner function Cw = Tr[/)Z}(A)], 
where p is the density matrix, and D{X) is the displacement operator. Accordingly the 
Wigner function is defined by 

TT^ J 

which at the origin then becomes 



ty(0,0) = At /rf^ATr 

TT^ J 



gAat-A*a^ 



Now using the Cambell-Baker-Hausdorf (CBH) theorem, e'^'^^ = e 2^^^^^e^e^ (as long 
as A and B commute with their commutator), we rewrite 
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W{0,0) = \ [(f\e~-^\^\\e^^^e-^'^), 



keeping in mind that the trace over the density matrix tensored with an operator is the 
expectation value of that operator. Expanding out the exponentials obtains 

/•oo /•27r p-^I^P / (—\\\'^)Pn'\PnP 

= d\x\d9\x\'-A^(y: ^ '^ " ^ 

io io ' ' ' ' 7r2 \^^ {p\y 

In the second line we have switched to polar coordinates (A = |A|e*^) and rearranged 
the sums into terms where p = q and terms where p ^ q. Now note that 

/ dOe''''' = 0, 
Jo 

for all integers x. So all the terms in the second sum integrate to zero and we have, 
after some rearranging, 

][ ^oo p2n 



w{o,o) = —Y. / d\x\de\x 



p 



{p\y I • 

The integral over |A| is not trivial but it can be changed to a known form with some 
substitution, the integral over 9 is straightforward, leaving 






2 



o oo n / r^\p 

^ n=0 p=0 P- 

where in the second line the expectation value is taken to be explicitly in the number 
basis (note that the sum over p can now only extend to n). The operators acting on 
the right number state produce 



^ „=0 p=0 \ P / 
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where all the factorials have been written as n-choose-p. Also note that we have 
multiplied by one. The sum over p is now clearly the binomial expansion of (1 — 2)". 
And so all that remains is 



o oo o 

w{o,o) = -J2Cn{-ir = -({-if 

^ n=0 ^ ^ 

The expectation value of the parity operator times a constant. Q.E.D. 
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